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ABSTRACT 


Cool star winds can best be observed in resonance absorption lines seen in the 
spectrum of a hot companion, doe to the wind passing in front of the blue stai. 

We calculate absorption line profiles that will be seen in the ultraviolet part of the 
blue companion spectrum. Line profiles are derived for different radial dependences 
of the cool star wind and for different orbital phases of the binary. 

Hovven and Willson (l!)S(>) find theoretically that stellar pulsations drive mass 
loss. We therefore apply our calculations to the Ceplieid binary S Muscae which has 
a B5V companion. 

We find an upper limit for the Ceplieid mass loss of M < 7 x 10 per year 
provided that the stellar wind of the companion does not influence the Ceplieid wind 
at large distances. 


1. INTRODUCTION 


1.1. The Problem 

In their studies of long period variables Bowen and Willson (198G) find large mass 
loss for these stars due to the pulsation, which expels mass during the expansion 
phases. Bowen and Willson conclude tentatively that all pulsating stars have strong 
mass loss anti that perhaps all mass loss of luminous stars is due to radial or noil-radial 
pulsations. 

I'or hot stars the mass loss can be best observed in the P Cygni profiles of the UV 
resonance lines of CIV at 1550/1, SilV at IWttA and MO.'M. For the cool luminous 
stars mass loss so far has mainly been studied in the Call K2 emission cores which 
may show blueshifted absorption components if mass loss occurs. For cool supergiants 
the Mgll resonance lines at 2S00/1 may also show absorption components in their 
emission cores (Stencel et. al. 1980) indicative of mass loss. For Gcphcids emission 
cores are rarely seen in these resonance lines and if they show absorption components 
we cannot be sure whether this means mass loss or whether the structure is due ot 
the pulsational motions. The observational determination of the mass loss rate is 
therefore difficult for these pulsating stars. 

Another method to determine mass loss of pulsating stars is provided by the 
observation of Cepheids with hot companions. We can then, in the spectrum of the 
companion observe the UV resonance absorption lines caused by the Cephcid wind 
passing in front of the* hot companion, see Figure I. Close to the companion the 
interaction between a possible wind of the companion and the Cephcid wind may 
cause some problems for the interpretation, these complications of the line profile 
occur near zero radial velocity components of the passing wind if observations are 



made at the phase shown in Figure 1(a). Large line of sight velocity components of 
the wind are only seen at large distances from the companion. Wind absorption seen 
at large blueshifts should not be influenced by the companion wind unless the latter 
is much stronger than the Cephcid wind. In that case we should see P Cygni line 
profiles of the companion in the observed spectrum. 

In order to determine the mass loss from the Cepheid binary S Muscae we will 
here calculate the absorption line profiles expected to be seen in the spectrum of the 
companion S Mus B, which is a B5 main sequence star. 

1 . 2 . Previous Tltcorcl.icnlStudics 

Absorption line profiles expected to be seen in the companion spectrum for hot 
star/cool supergiant binaries due to the suporgiant wind were previously calculated 
by Ilempe ( 1 OS'I ) with the o Sco system in mind. lie also included the emission by 
the circumstellar wind material in his line profile calculations. We calculate emission 
for one set of parameters only and show that for our application it is not important. 
Ilempe considered a velocity law 

v{r) = Voo{l - II a/ r) 05 , (I) 

where Ra is the radius of the mass-losing star and Voo is the wind velocity at large 
distances from the star A. Ilempe claims that if the turbulent velocity of the wind 
material is comparable to the wind velocity the exact form of the velocity law does not 
matter unless the hot companion is ‘behind’ the cool star. Under most circumstances 
all velocity laws will produce very similar absorption line profiles. 

Our calculations show that this does not hold for velocity laws that decrease 
with increasing distance from the mass losing star and also velocity laws which give 
a maximum velocity at some point. While Ilempe explores a large range of optical 
depths along the line of sight and also a large range in r>oo/ l ’turl> W( - restrict our 
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calculations to small optical depths and only a few values of VtmV (Here, u tur b is the 
turbulent velocity in the wind.) Instead, we concentrate on the changes in absorption 
line profiles due to different velocity laws. 

The organization of this paper is as follows. In Chapter 2 we go over our methods 
of calculating line profiles. In Chapter 3 we describe the absorption lines calculated 
for different velocity laws and different orbital phases. We also show the emission 
line obtained for one special case. In Chapter 4 we discuss the observations for 
the S Muscae system and compare theoretical and observed line profiles in order to 
determine the mass loss for S Muscae. Chapter 5 contains a. discussion of errors and 
finally Chapter (i summarizes our results. 


2. METHOD OF CALCULATION 

2.1. The Geometry 

The geometry of the binary system is shown in Figure 1. The coordinate parallel 
to the line of sight is z, and 2 = 0 in a plane perpendicular to the line of sight that 
contains the Cepheid. The two stars describe circular orbits around their common 
center of mass; the distance between the two stars is assumed to be constant. Figure 
1 also explains our notation for the orbital phases, <j). Angle (f> is the angle between 
the line of sight and the line going from the Cepheid to the L3 star. 

It can be seen from Figure I that absorption lines sample diflerent columns 
through the wind at different orbital phases. 

For position depicted in Figure 1(a) only material with radial line of sight, velocity 
u r < 0 is seen. The wind is sampled for distances r from the cool star for which r > d 
where d is the orbital radius of the blue star around the cool supergiant. 


For position in Figure 1(b) Lite absorption lines sample the wind for distances 
r > 6 where 8 < d. Distances r < d from the cool star are seen twice: once with the 
positive and once with negative line of sight velocity component v r . 

In position shown in Figure 1(c), when the hot companion is ‘in front’ of the cool 
star, only distances r > d are sampled by the absorption line. 


2.2. Calculations of Absorption Component 

In calculating model line profiles we make some simplifying assumptions about the 

nature of the wind similar to the ones made by Ilempe (1984). Firstly, mass outflow 

from the mass losing star A is assumed to be spherically symmetric about the center 

or the star. In a real situation there are a number of factors that can introduce 

asymmetry into the wind, for example, gravitational perturbation of (.he mass flow 

by the hot companion star B, radiation pressure from star B, clumps resulting from 

instabilities causing small local shocks and turbulence, rotation of the mass losing 

star, orbital motion of the binary and wind from star B (if present). We do not 

include these complications. Velocity and density profile as a function of distance 

from the star are assumed to be constant in time, this means, for instance, that 

pulsations of the Cepheid are ignored. This may be an unrealistic approximation for 

this special system since the pulsations may be the very mechanism responsible for 

driving mass loss. We will come back to this question in Chapter V. 

_ v.fU Y iX of \‘A\ y 'k' 

As will be explained later onIy\high T velocity components or <lhe 

pcVvio , ' J 

studied in detail and modeled. These wind components arise far from the two stars 
of the binary and hence we feel justified in assuming that ionization fractionjof ions 1 

like Fell, Sill, A II II, etc. remaiiuf constant. j 

The absorption along the column ‘in front’ of star B (see Fig. I) is calculated as 
follows. The column is divided into a large number of segments of length dz. The 
optical depth d.T\ is evaluated for each segment, and then contributions due to all line 


r ,uU ' pU lc 
wind^ will be 
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segments dz are added together to get the total optical depth of the line of sight at 
one particular wavelength, 


dr x = n(r) ■ k q • c" (AA / AA,,)2 • dz 


0 >) 


Here, r is the distance of the line segment dz from the mass-losing star A, »(r) 
is the number density of absorbing particles in the line segment dz, and /<o is the 
absorption coefficient per absorbing particle in the center of the line in the rest frame 
of the absorbers, and 


A A o = A n 


^lurb 

r 


( 3 ) 


Turbulent velocity, n Jin | M is n. parameter to be determined by comparison with the 
observations. The length of each line segment dz was chosen such that the change in 
the line of sight velocity across dz is very small compared to the turbulent velocity 
in the wind. 


Conservation of mass in the wind requires that 


47rr 2 N(r)v(r) = constant. 


0 ) 


For a given v(r) this determines total (i.e., all species) number density N(r). The 
abundance and the degree of ionization of the line producing element determines 
n(r)/N(r). 


2.3. Calculation of Emission Component 
The bulk of the emission comes from a large region (linear size of about twice 
the stellar separation) around the binary. Emission from this circumstellar envelope 
is dominated by radiative excitations from the light of the B star and deexcitations 
- collisions can be shown to be unimportant. 

We approximate the radiation field of the companion by a Planck function for the 
T c (r of the hot companion B. The number of excitations per unit volume and second 
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from the ground level at distance r# from star B is then 


^exc — Ji 


( 5 ) 


where 


J v = \b v (T^)^. 


4 ' vq 1 

Bi u is the Einstein B coefficient for stim ulated excitation from lower level / to the 
upper level u. The energy emitted into a unit solid angle in the line is given by 


e = N exc Iw/4n. 


(G) 


The wavelength dependence of the emission in the rest frame of the emitting part icles 
is determined by Doppler shifts due to turbulent motions and is given by 

V?(AA) = (7) 

and is centered at A r = A«(l -f ^), if v r is the line of sight velocity component of the 
volume element under consideration. 

The energy emitted is reduced due to absorption along the path through the 
wind. This absorption is calculated in the same way as the absorption in front of 
star B. Integration of the emission c at each wavelength AA = A — A () over all volume 
elements divided by D 2 gives the amount of energy received here per cm 2 per second. 
D is the distance between the Cepheid and us. 

The energy which we receive from star B per cm 2 and second is given by 

F x = IWd r)^r«- n - ( 8 ) 

Both fluxes are proportional to (Rd/D) 2 . In the normalized profiles these factors 
cancel. 

The emission line profiles calculated here look similar (in general) to ITcmpe’s but 
are not completely identical to his profiles. 
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2.4. Calcula tion of Mass Loss Rate 
The mass loss rate from the surface of star A is 

M = At: • No • 7 'o • m//. (9) 

v 0 is the wind velocity at the base of star A and mjj is the mass of Hydrogen atom. We 
assumed that metals do not contribute much to the mass of the outflowing matciial. 


3. ABSORPTION LINE PROFILES FOll DIFFERENT VELOCITY LAWS 

3.1. The Standard Velocity Law 

The velocity law most commonly adopted for luminous OB stars, A-type super- 
giants (Castor k Earners 1979), and stars with dense expanding winds, like Of and 
Wolf-Rayct stars (Olson 1981) is of the form, 

v(r) = i>oo(0.01 + 0.99 [1 - ly-f). ( 10) 

Exponent f3 is usually assumed to be of the order of 1. Figure 2 shows the radial 
dependence of the velocity according to equation (10) with /? = 0.5. 

The corresponding absorption profiles for phase angles <j) = 30°, 60°, 90° and 150 
(see Fig. 1) are shown in Figure 3. The scale on the horizontal axis is in terms of the 
terminal wind velocity, v (X . Ratio Woo/^turb = 10. I he corresponding mass loss rate 
if the line is A23S2.03 of Fell is approximately M ~ 1.12 x lO“ lo A/ 0 .yr -1 . 

'I’he absorption at the shortest wavelengths corresponding to the largest radial 
velocities arc due to the material at the largest distances from star B. Looking at. 
the absorption line profile in Figure 3 it seems at first sight surprising that the line 
depth increases for larger radial velocities v T even though the wind density decreases 
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for larger radial velocities i.e. larger distances from star A. It can, however, be easily 
verified that the increase in column length contributing to absorption between v r and 
v r + dv r more than compensates for the decreasing density. 

The sharp dip at the shortest wavelengths is due to material at terminal velocity. 
Even though the density of absorbers with large u r , i.e. at large distances from the 
star is small, the length of the column contributing to absorption is theoretically 
infinite. 

The sharp dip at terminal velocity is common to all the profiles because the 
absorption at large distances is not affected by the phase <J) of the system. For phases 
<f> > 90° there is no absorption on the long wavelength side of Au because star B is ‘in 
front’ of star A and only the part of the wind which is approaching us is seen in front 
of B. Phases with small <jS’s (for example, <j> = 30°) give rise to deeper absorption 
profiles than for larger <f> because at small 4> our line of sight to star B passes close to 
star A where the density is high and the wind velocity small. Comparing absorption 
profiles at 4> and 180° — 4>, for example, </> = 30° and 150°, one can estimate the amount 
of absorption in a column between the two positions of star B. This information is 
useful in mapping out the spatial and velocity structure of the wind. 

The ratio of the turbulent velocity to the terminal velocity of the wind has some ef- 
fect on the shape of the line profile. Turbulence naturally broadens the line, smoothens 
out sharp edges and spikes. It also moves the position of the lowest flux in the absorp- 
tion profile (like the one in Fig. 3) to longer wavelengths. This happens because some 
of the absorbers moving at terminal velocity will be ‘redshifted’ by large turbulent 
velocity thus making the number of absorbers at velocities just shy of the terminal 
velocity large. The solid line in Figure 4 is the same as the <f> = 90° profile of Figure 
3. The other two lines of Figure 4 difier from the solid line oidy in the ratio Uoo/ ?, i,i r l> 
as shown. 
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3.2. Exponential Velocity Law 

In the last section vve showed that a monotonically increasing velocity law gives 
rise to a profile with a broad shallow shoulder and a deep spike-like part on the short 
wavelength side of Ao- This spike has an abrupt blue edge. How much and in which 
way does the absorption line profile change for different velocity laws? 

In order to study this we modified the wind velocity law such that after an initially 
high velocity close to the mass losing star the velocity drops at large distances from 
the star. This is achieved with an exponential velocity law of the form 

v(r) = eo(<t • r~’^ -f r) (II) 

A typical example of such a velocity law is shown in Figure , r >, with « — Iff), b = 2.5, 
c = 0.12 and e 0 = 100. Note that in this particular case the velocity nearly reaches 
its asymptotic value at r/ H a ~ 15, which is larger than the separation of the stars 
in the S Mus system. 

Since mass flow and radiation field come from two spatially separated sources, 
the radial velocity distribution, v r , of the matter in a column in front of star B (i.e., 
material responsible for absorption) can be quite different from n(r). This point is 
well illustrated in the present case of an exponentially decreasing velocity law. Figure 
6 shows the radial velocity profile along the line of sight as a function of 2 . There is 
a hump in v r , even though v(r) is decreasing monotonically. 

Absorption line profiles for this case arc shown in Figure 7 for different orbital 
phases. The velocity is expressed in terms of the terminal velocity, which is v^c (see 
eq. [11]). For small angles <j> our line of sight to star B passes close to star A where 
the velocity of outflowing material is large, so the red and blue wings on either side of 
A 0 are farther away from the rest wavelength than for larger values of <f>. Conversely, 
for 4> close to 90°, the separation of the two stars on the sky is large. Now light from 
star B on its way to us passes through slow moving material at large distances from 
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star A and hence all absorption is very close to A 0 . The central dip which arises due 
to the material at terminal velocity, is not alTected by the phase of the system. 

We discuss in detail the absorption line profile generated by this velocity law 
at orbital phase <f> = 30°. The dip just on the short wavelength side of the rest 
wavelength is due to the wind material at large distances from the star. In addition 
to this main feature there arc also two less prominent features — - the shallow broad 
wings on either side of the line center. The wing at the blueshifted wavelengths comes 
from material at and near the line of sight velocity peak. Even though the density is 
quite low there because of the high velocity and large distance from star A, there arc 
two spatially separated regions in the line of sight contributing to the absorption at 
those e r ’s (see Fig. (>). A similar wing appears on the red side it is due to material 
from the ‘back 7 side of star A which is streaming away from us. 

The main difference of this type of absorption profile compared to that of equation 
(10) (see Fig. 3) is that the position of the maximum absorption dip is changed. In the 
case of the standard velocity law most absorption occurs at the shortest wavelengths 
while in the case of the exponential law it is very close to the rest wavelength. In 
both cases the amount of absorption at velocities other than terminal velocities varies 
with the orbital phase of the system. 

3.3. Velocity Laws with a Maximum 

This is a modification of the exponentially declining law considered in the previous 
section. Equation (11) implies that velocity of the material at the surface of star A 
is very large. For example, taking equation (11) with a = ill), b = 2 . 3 , c = 0.12 
and = 100 gives 2023 kin.s -1 as the velocity close to the surface of star A, where 
r = R,\. This rather large velocity is very difficult to explain in physical terms since 
it means that the wind is highly supersonic at the surface of a star. A more physically 
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plausible situation is achieved by combining equations (10) and (11), namely 

v{r) = «n(0.01 +0.99(1 - ^]) x (a • e r r ! R * h + c). (12) 

It can be easily seen from equation (12) that now the initial velocity is smaller by a 
factor of 100 (because 1-^4 is small there), i.e. v(Ra) is now ~ 20 krn.s” 1 . The 
velocity then rises rapidly, reaches a peak close to star A (at vj Ra ~ 2) and declines 
exponentially as in equation (11), see Figure 8. The velocity according to equation 
(10) dominates the one of equation (12) only at small r. The velocity according 
to equation (12) is therefore similar to the velocity of equation (11) at most radii 
and therefore the resultant absorption lino looks very similar to that obtained with 
equation (II). (see Fig. 0.) Note, however, that the densities at the surlace of the 
star, and hence the mass loss rates in tin* two profile's of Figure* 9 are very different! 

This wind profile is qualitatively similar to the wind profile of some cools stars 
as given by Carpenter (1984). lie deduced the shape of the velocity profile from 
the Fell asymmetric line profiles in the extended chromosphere of a Orionis. I he 
strongest parts of the line are formed farthest from the photosphere in a decelerating 
region, intermediate strength line parts are formed in accelerating regions closer to 
the photosphere and the velocity in the part of the wind closest to the star is almost 
flat. 
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4. THE S MUSCAE SYSTEM 


4.1. The Geometry of the S Muse ae System 

The physical parameters of the S Muscac system arc as follows. S Mus is a binary 
with an orbital period of 505.2 days and a mass sum of approximately 10M© (Evans 
k Bolton 1990). The pulsational period of the Cepheid is 9.66 days (Lloyd Evans 
1982). With an observed orbital velocity amplitude u r sim w 13.8 km.s -1 we derive 
an inclination i « 30°. Star A, which is a Cepheid, has a radius of roughly 50/2©. 
Radius of star B, a B5V star, is approximately 5/?©. Our assumption of a circular 
orbit for the stars in the binary is justified here since the eccentricity ol tin* S Mus 
system is, e = 0.03 db 0.03 (Lloyd Evans 1982, Evans k Bolton 1990). The two stars 
are separated by ~ 576/2©. The effective temperature of star B was taken to be 
17, 000° K (Bolim- Vi tense k Proffitt 1985), and its spectrum in the UV region was 
approximated by a black body at that temperature. 

The geometry of the system and the direction to the observer at different orbital 
phases arc shown in Figure 1. As was explained in Chapter 2, at different orbital 
phases the line of sight to the hot companion may pass through the wind at different 
distances from the mass losing star. 

In addition to determining the wind velocities at different distances 7- from the 
cool star we can also sample the degree of ionization of the wind at different distances 
from the star by observing the line profiles for lines of different ions. Such studies are 
only possible in the ultraviolet where strong resonance lines of Eel, Fell, Mull, Sil, 
Sill, Sil 1 1, SilV, Cl I, CIV, All I and A 1 1 1 1 can be observed in the spectrum of the B5 
star unimpeded by the varying Cepheid lines. Note that the degree of ionization of 
the cool star wind may be expected to be higher close to the blue companion either 
due to a shock developing around the B star if it accretes material or, if the B star is 
losing mass, due to the shock in the region where the winds from both stars collide. 


(See Che-Bohnenstengel & Rciiners 198G.) 

4.2. The Ultraviolet Spectra of S Muscae 
The observational data used in this paper was taken by the IUE satellite in the 
period 1982-1988 using LWR, LWP and SWP cameras. Most of the observations were 
done by one of us (E. B.-V.). 1 he complete list of the spectra, the dates they weie 
obtained, the orbital phase of the S Mus binary system and the pulsational phase of 
the Cepheid at those times are given in Table 1. Most of the spectra were taken at 
the maximum radial velocity phase, i.e. orbital phases 0.25 and 0.75. 

Figures 10, II, and 12 show spectra of several ions at various orbital phases of the 
system (see Table I). (Note: only standard IUE processing was used to obtain these 
spectra.) All spectra in Fig. 10 and I I were taken at an orbital phase <f> = 90°. The. 
line profiles for A < 2400/1 are very noisy because of the low sensitivity ol the camera 
for these wavelengths. All the resonance lines have a very deep absorption feature 
dose to the 0 radial velocity. 'Phis central feature does not, vary with the phase of the 

iVn>(mUt«. e. (a/Y^-C , . 

system and -fs believed-te-be mostly due to the interstellar material V the equivalent 

widths of these central lines are somewhat larger than what would be expected given 

the reddening of S Mus. The widths of these lines (^FWIIM ~ 28 km.s 1 ) arc 

-V>uU Ciclcv-CA cooTaIU'-Iauv’ f<tW yUw, • , 

slightly higher than the typical velocity dispersions of interstellar clouds! IP part of 

the central feature is due to the wind of star B then that could in principle account 

for the width of the line and also the fact that <f> - 90° lines are slightly redshifted ’ . 

my-'''' ;,f 

from 0 km.s 1 . « r , . ... TA 

At orbital phase <f) = 90° most of the 0-volt resonance lines of Fell (Fig. 10), Sill, 

A1I1I with moderately good signal to noise ratio, have a gentle extended dip at about 
-170 km.s -1 from the line center. We need a theoretical profile that reproduces this 
blue dip at <f> = 90° and possibly shows some absorption close to the rest wavelength. 

The extended blucshiftcd dip is absent in some other resonance lines, most notably 
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in the weaker Mull lines (Fig. 11, 12). The absence of this feature in non-resonance 
lines can be explained by the lower abundance of Mn as compared to Fe, about a 
factor of 100. 

There is another absorption feature which is common to many lines. The short 
wavelength side of the central absorption dip has a gentler slope than its long wave- 
length counterpart. This feature is quite prominent in the lines of Fell, Mull, SilV, 
A1III at all orbital phases of the system (Fig. 10, 11, and 12) though it appears to 
be variable. 

This feature shows low line of sight velocity components and therefore could pos- 
sibly be an interstellar component. In order to check this we looked at the spectrum of 
a 132V star in a different direction of the sky, IID 3309, and found similar line profiles 
for several lines. This suggests that this extended blue wing is probably intrinsic to 
the star. The apparent variability of the feature argues against, an interstellar origin. 

This feature might then be due to a 13 star wind. The appearance of this extended 
wing also in the supposedly single star IID 3369 would argue in favor of this inter- 
pretation. We checked the line profile which would be expected for a single B star 
wind assuming the velocity law for this wind to be of the form given in equation (10). 
Figure 13 shows the expected line profile, which looks rather similar to the observed 
one. (Note, however, that in this case the total number density at the surface of the 
single star has been reduced by a factor of 10 to make the profiles look similar.) It 
surprises us, however, that the wind of a star with 17,000K should have such a low 
degree of ionization. We do however also see some indication of a ‘blue’ wing in some 
higher ionization particles like CMV. 

It may also be possible that this feature could be due to the Cepheid wind close 
to the B star. In any case the interpretation of this feature is ambiguous. We there- 
fore concentrate on the shorter wavelength absorption components which originate 
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in material with a large line of sight velocity which means material which is at large 
distances from star B. This material is not disturbed by the B star wind and gives 
us more reliable information about the Ccpheid wind. (If the B star wind carries a 
larger momentum than the Cepheid wind it could blow away the Ceplioid wind, we 

' |M l 

would then see a typical stellar wind line profile for the B star lines, which we do not 
observe.) 


4.3. Fitting Model Line Profiles To Observations 

We have the best spectra and also most spectra for phase (f> = 90°, We shall 
therefore use these spectra. The shallow dip on the blue side ol the rest wavelength can 
be well approximated by the monotonically increasing velocity law given by equation 

V/ N 1 » '»■**'* ■ 1 S ft • b : f ( v I t \ : •. ’e ' ' - . 

(10). The other two velocity laws’ would yield extended dips' (‘satellite’ features) only 

'V 1 '•/' » 

at phases other than = 90°. (See Fig. 7, 9.) Therefore we arc going to assume the 
velocity law of equation (10) with /? = 0.5 in all our subsequent profile comparisons. 

Theoretical profiles based on this velocity law and the geometry of the S Muscae 
system were fitted to several Fell, Mull and A 1 1 II lines. (Fig. Tt) Note that it is only 

the shallow blue wing around -170 km.s -1 that is being modeled here. The central 

{ / 

C v • • * V*A,f e { , ;■ * M < r • r. ■ I • • l . » 

feature (dose to 0 km.s -1 ) is probably mostly interstellar as discussed earlier. Using 
solar abundances of elements, i.c. Fe/II=3xl0 — 5 , Mn/II=2.5xlO 7 , Al/II=2xlO G , 
the set of parameters that fit the observed profiles best are as follows (solid line): ?>«, = 
220 km.s -1 , u tll rl) = 50 km.s -1 , and the mass loss rate, M = 1.5 x 10 -lo M©yr -1 
assuming that Fe is all in the form of Fell. Since the ‘dips’ differ a. little from line 
to line and from spectrum to spectrum due to mainly noise, there is a range of 
parameters that can fit the observed line profiles equally well. For example, the 
following set of parameters will also fit the profiles (dashed line in Fig. Tl) t?oo = 185 
km.s -1 , u tHrb = 37 km.s -1 , and M = 1.13 x 10 -,0 M©yr -1 . (Note that lines whose 
spectra do not show any blue wing absorption (e.g. A IT 1 1 AIS5T72, All II A 1862.78 
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and Mnll A2605.70) are not expected to have any observable absorption according to 
the model.) The mass loss rates in both cases are very similar. The mass loss rate is 
Unis smaller than might have been expected for a pulsating star if mass loss is indeed 
mainly occuring due to pulsations. Could the emission from the envelope influence 
the result? We have calculated the emission from the envelope around the binary 
using the second set of parameters above. Figiuc 15 shows a 2382/1 resonance line of 
Fell with four theoretical lines superimposed on it. A pure absorption line and a final 
line profile including emission. The two types of line, solid and dashed, have the same 
meaning as in Figure 14. A slight improvement in the fit may perhaps be seen when 
including the weak emission. The emission part of the line profile is not, however, 
important in comparison with the absorption for the following reasons: firstly, the 
strength of the emission was always found to be much less than the corresponding 
absorption, and secondly, since the bulk of the emission comes from the region between 
the two stars where both the density of the wind material and the radiation intensity 
are high, the emission is usually not far displaced from the line center where it is 
almost completely covered up by much stronger interstellar absorption. Because of 
these factors we have concentrated on the far wings of the absorption corresponding 
to high radial velocities. 


5. ERROR DISCUSSION 

Several sources of error have to be considered. 

5.1. Placement of the Continuum 

Since the short wavelength absorption wing is wide ami shallow and since the 
spectra are noisy we may misjudge the position of the continuum. Because these 
wings are optically thin the derived mass loss rate is directly proportional to the 
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depth. The true line depth would have to be twice as deep as adopted for the mass 
loss rate to be twice the derived value. This seems to be an upper limit for the 
uncertainty in the line depth. An error of this kind will not change the ordei of 
magnitude for the mass loss rate. 

Any error in the adopted oscillator strength would also be of minor importance. 

5.2. Degree of Ionization of the Wind 

In deriving the mass loss rate we assumed that iron is all in the form of Fell. 
There are no usable Fell! lines in the II JR wavelength region. We therefore looked 
at the AIIII resonance lines at 1854.7/1 and ISG2.S/1 to see whether any blue wings 
due to mass loss could be seen. Figure M shows the observed line profiles. We see 
no indication of blue wings. All I has an ionization energy of 18.S2cV while hell has 
an ionization energy of IG.lSeV. The difference is very small. The abundance of A1 is 
a factor of 10 less than Fe. The oscillator strengths of the AIIII lines are a factor of 
2.G larger than the one for the 2599.4/1 line (/= 0.5G5 for the AIIII 1854.7/1 line and 
/=0.22 for the Fell 2599.4/1 line). If there would be 5 times as much AIIII as All I 
we should see the shallow blue wing of AIIII line as strong as the Fell line. Since we 
do not see any AIIII we conclude that there could be at most 3 times as much belli 
as Fell if any at all. 

We thus conclude that the derived mass loss rate due to these effects could be 
wrong by at most a factor of G, if both the uncertainties conspire to work in the same 
way. 


5.3. Reality of Blue Wing 

We might ask whether the shallow dip around -170 km.s" 1 is actually due to the 
wind of the companion. Could it be intrinsic to the B5V star spectrum? The profile 
does not support this suspicion, however we did look at the UV spectra of IID G0098, 
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a standard single B4V star. Sonic lines in its spectra seemed to liave a shallow dip at 
around -170 km.s -1 . There are very few good IUE comparison spectra for standard 
stars of this spectral type. Most of the observed stars are either peculiar or have high 
rotational velocities. In any case assuming that the shallow dip is due to the Ccpheid 
wind gives us an upper limit for the Ccpheid mass loss. We then conclude that for 
the S Mus Cepheid, M < 7 x 10 _,o A/© per year. 

5.4. Periodic Variations 

Could it be that the wind of the Cepheid varies periodically and we just happened 
to study the wind absorption at a phase when little material is passing by the B5V 
star?' For the phase </> = 90° wo have observations for a whole pulsational cycle (Bdhin- / 

Vitense el. al. 1990). We do not sec' any systematic variations for this time interval. 
Should we expect any variations at all? 

If we assume that the Ccpheid mass loss occurs in discrete bursts then we might 
be able to see these in the absorption line. Along the line of sight we should see density 
enhancements at certain intervals. These high density regions would correspond to 
different line of sight velocity components. We might then sec increased absorptions 
at regular intervals. We find, however, that the line of sight velocity separations 
of these high density regions are typically smaller than the turbulent velocity in the 
wind. Hence any discreteness in the wind that might have been there to start with will 
be washed out by the turbulence in the wind. So we do not expect to see variations 
in the shape of the absorption line during the course of a pulsational cycle. 

5.5. Disturbance of the Cepheid Wind by a Possible Wind of the BhV Companion 

Little is known about winds from B5V stars, but they are assumed to be weak and 
the corresponding mass loss small. (Note, however, that Hagen, Ilompe and Renners 
(1987) claim to see a wind from a Sco companion (B4V) with about 500 kin.s -1 in 
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the Sill lines, which is however not seen in the Fell, CrII or Nill lines.) 

If a B5V stellar wind is present its mass loss line profile should not vary in any 
systematic way for different orbital phases while the Cepheid wind profile should vary 
systematically in the manner described above. The two stellar winds should thus be 
distinguishable. 

For the S Mus companion we do sometimes see an extended wing in the SilV 
1393A line, but it is not seen in the weaker 1402/1 line. It is therefore likely that the 
wing in the 1393/1 line is due to another line. We may see an asymmetry in the A 1 1 1 1 
1854.7/1 line and in the Fell line at 2585,9/1 (Fig. 14) which may be due to a B5V 
star breeze. 

The asymmetry in the line extends to e r ~ —50 km.s" 1 . For a B5V star the 
escape velocity is ~700 kirns"' 1 . It is therefore doubtful whether the asymmetry is 
indeed due to a wind which escapes from the star. If it is due to an escaping wind 
then the high velocity densities are so small that we do not see most of the wind 
absorption line, we only see the part of the line at small velocities, v T > —50 km.s -1 
, i.c. dose to the B star. This means that at velocities v r ~ 1 00 km.s' 1 which we 
study (p • u) from the B5V star must be less than (p • v) from the Cepheid, where p 
and v are density and velocity of the wind respectively. The momentum from the B 
star wind is therefore smaller than the one from the Cepheid. We therefore do not 
think that it is possible for the B5V star wind to "blow” the Cepheid wind out of the 
line of sight column. 


6* CONCLUSIONS 

In this paper wc have considered three different types of velocity laws for mass 
loss from stars. The three types arc: a monolonically increasing law, a monoloni- 
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cally decreasing law (exponential decline), and a combination of both, a wind which 
accelerates, reaches a peak and decelerates to a constant terminal velocity. 

The corresponding absorption lines arising in the spectrum of the companion were 
calculated for different parameters of the velocity law and different orbital phases of 
the binary system. The material responsible for absorption on the one side and the 
radiation field on the other side were assumed to come from two spatially separated 
sources, star A and star B respectively. This special geometry gave rise to interesting 
absorption line shapes. 

l-'or monolonieally increasing wind velocities we find the deepest absorption at 
the short wavelength edge of the absorption line at the terminal velocity. 

| '’or outward decreasing wind velocities the radial velocity along the line ol sight, 
shows maxima leading to more or less broad and fiat satellite components whose 
wavelengths depend on both, orbital phase and the velocity law. Observing the 
absorption lines at different phases permits the separation of the phase effects and 

the wind velocity law. 

The analysis of the spectra of S Mas system is complicated by a number of factors: 

, t , Ui, i iv c* ^-«d> v , 

firstly part of the observed absorption featured arcTdue to interstellar material, and 

liU/ l W ^ tlri Gl *■ d l i* \a •£ s tT l it* t c iv i r\ b tt i J 

part could be due to the wind of the hot companion*. Secondly, the gravity field of stai 
B could, in principle, influence the pattern of the wind material velocity, such that 
part of the material close to the B star might be falling towards star B. In this case we 
should have some absorption on the long wavelength side. There is some indication 
of that in the spectra of some lines but the evidence is insufficient. Thirdly, most of 
the spectra are plagued with irregular bumps and dips, some of which arc, ol course, 
noise, but some could be due to the non-uniform nature of the wind. For example, 
turbulence in the wind can produce local shocks which will lead to a distoited line 
profile shape. By observing only the parts of the wing for v r < — ' >0 kirns" 1 we avoid 
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problems with material close to the 135 companion. This also avoids the problems 
with the interstellar lines. 

In order to deal with the above difficulties, we have to make some simplifying 
assumptions about the S Mus system and then proceed with the fitting. We assume 
a spherically symmetry wind which should be a good approximation for the large 
distances from the B5 star which we are studying. The velocity law that fits the 
observations best is given by equation (10). The corresponding mass loss is derived 
to be of the order of 10~ 9 - lO _lo A/0yr _1 . This would be an upper limit if the shallow 
blue wing dip were not entirely due to the Cepheid wind. Considering uncertainties 
in the continuum lil ting and in the degree of ionization the mass loss could possibly 
bo larger by a factor ol G. 

We do not think that a possible wind from I35V star companion could disturb 
the Cepheid wind enough to have an effect on our mass loss determination. 

The main conclusion of this study is then that whatever the exact form of the 
velocity law is and whether or not star 13 has a wind of its own, the mass loss rate from 
the Cepheid is very small, < 7 x lO -1( M/ ( r ) yr - ' 1 , A more comprehensive study is nec- 
essary to determine whether mass loss from other, especially longer period Cepheids 
may be larger. A search for long period Cepheids with suitable blue companions is 
necessary to make such a study possible. 
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TABLE l 


Observational Data 


Spectrum 

Date 

Orbital Phase 

Pulsational Phase 

LWP 3035 

1984, 7S day 

0.29 

0.030 

LWR 13849 

1982, 215 day 

0.10 

0.700 

LWR 15538 

1883, 81 day 

0.56 

0.625 

LWR 15975 

1983, 139 day 

0.67 

0.643 

SWP 31542 

>> 

•M 

CM 

r- 

00 

O 

0.72 

0.739 

SWP 31543 

1987, 224 day 

0.72 

0.739 

SWP 33435 

1988, 124 day 

0.25 

0.125 

LWP 13163 

1988, 124 day 

0.25 

0.125 

LWP 13176 

1988, 12S day 

0.25 

0.562 

LWP 13183 

1988, 129 day 

0.25 

0.643 

LWP 13196 

1988, 131 day 

0.25 

0.944 

LWP 13233 

1988, 136 day 

0.25 

0.436 
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Figure Captions 


Fig. 1. — The geometry of the binary system at various phase angles, <j>. A 
is the mass losing star, and B is the blue companion. The shaded column shows 
the material seen in absorption by the observer, (a) <^=90°, (b) <fi <90°, (c) 
(f> >90°. 


Fig. 2. — Standard velocity law as a function of the distance from the mass 
losing star A. The form o f the velocity law is v(r) = (0.01 + 0.99[1 — Ra/ 1 '] 0 ' 5 )- 

Fig. 3. — Absorption line profiles from a binary system. The velocity law is 
shown in Fig. 2. Different line types correspond to different phase angles of the 
system: </>=30°, 60°, 90°, and 150°. The ratio of the terminal wind velocity to 
the turbulent velocity in the wind, i>oo/ J ’turln i- s equal to 10 for all four profiles. 

Fig. 4. — Absorption line profiles for the velocity law shown in Fig. 3. The 
three different types of lines correspond to three ratios of the terminal velocity 
of the wind to its turbulent velocity: t’oo/ftiiib—T 10, and 15. 


Fig. 5. — Exponential velocity law as a (unction ol the distance from the 
mass losing star. The form of the velocity law is v(r) = vo(30e -r / 2,5/</ * + 0.12) 
km.s -1 . 


Fig. 0. — Radial profile, t> r , of the exponential velocity law shown in Fig. 
5 as a function of the radial coordinate, z. The phase of the binary is <^>=90°. 
Notice that the profile is not a monotonous function; in particular, it exhibits 
a peak at a distance of 5/?,i from the blu e com panion. 

Fig. 7. — Absorption line profiles at different orbital phases of the binary, 
<^=30°, 40°, C0°, 90°, and 150°. Velocity law shown in Fig. 5 was used. The 
ratio Voo/i’uirb > s 10 for all five profiles. 

Fig. 8. A velocity law with a maximum as a function of the distance from 
the mass losing star, u(7 - ) = f>o(0.01 + 0.99[1 — Il/\/r]) • (30e -r / 25 ^ /1 + 0.12) 
km.s -1 . Notice that the velocity profile has a peak at a distance smaller than 
the separation of the two stars in the binary. 
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Fig. 9 # — A comparison of absorption lines profiles for two different velocity 
laws: the solid line corresponds to velocity law of Fig. 5, the dashed line 
corresponds to that of Fig. 8. The orbital phase (f> is 30° in both cases. Notice 
that even though the shapes of the two line profiles are similar in both cases, 
the number of absorbers along the line of sight is not. Mass loss rate for the 
dashed line profile (velocity law with a maximum) is 2.5x10° times larger than 
the one for the solid line profile (exponential velocity law). 


Fig. 10. — Observed line profiles of five resonant Fell lines. All spectra were 
taken at an orbital phase <p=\)0° . 


Fig. 11. — Observed line profiles of four resonant Mull lines. All spectra 
were taken at an orbital phase </>=90°. Note that four features in this Figure 
are crossed out and labeled /?. These are reseau marks. 


Fig. 12. -Observed line profiles of resonant lines of Mull. These spectra 
were taken at orbital phases other than </>=!M)°. 

Fig. 13. — Absorption line profiles of the velocity law given in Fig. 2. The 
solid line shows the profile that would arise from a single star. 'The dashed line 
corresponds to a line seen in a binary, where mass outflowing from one star is 
illuminated by light Irotn the hot companion star. 

Fig. 14. — Six panels show observed line profiles of several resonant lines 
and two model fits to each line. Note that it is the shallow wing at around -170 
km/s that is being fit, not the central part of the line! Both models assume 
the standard form of the velocity law shown in Fig. 2. The parameters of the 
velocity law are as follows: solid line, Uqo = 220 km.s -1 , Vt U rb = 50 km.s -1 , and 
the mass loss rate, M = 1.5 x 10 -1(, iVT y r ^ ; dashed line, Uoo— 185 km.s 
*Aurb=37 km.s -1 , and M = 1.13 x 10 -10 jW- yr -1 . Note that we did not attempt 
to model the central strong component of the absorption line because it is 
probably mostly interstellar. 

Fig. 15. — The influence of the emission part of the line on the line profile 
as a whole. A resonant line of Fell is used as an example. One set of solid and 
dashed lines represents absorption only, while the other set of lines corresponds 
to the absorption plus emission. It can bo seen that emission introduces very 
little difference into the line profiles lar from the rest wavelength. (Note that it 
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is only the shallow blue wing that is being fit.) 
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